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I . INTRODUCTION 

In  previous  reports  the  details  of  the  Molecular  Stuffing 
process  have  been  discussed.  Briefly,  an  alkali  borosilicate 
glass  is  melted  by  conventional  means,  formed  into  rods,  phase 
separated  and  leached  to  produce  a porous  structure  which  is 
then  doped  with  an  index  modifier  and  sintered  to  yield  a pre- 
form. An  object  of  the  research  under  this  contract  is  to  de- 
velop the  ability  to  control  the  doping  to  produce  a nearly 
parabolic  index  profile  in  the  preform. 

This  final jeport  for  FY-1976  presents  the  results  of  dif- 
fusion coefficient  measurements  for  doping  solutions  in  porous 
rods  and  results  of  a computer  study  done  to  determine  values 

i 

for  the  parameters  necessary  to  produce  a parabolic  index  pro- 
file by  Molecular  Stuffing. 

II  PROGRESS  IN  THE  DEVELOPMENT  OF  A STUFFING  PROCEDURE  TO 

I 

PRODUCE  PARABOLIC  INDEX  PROFILES 

If  one  knows  the  composition  profile  (corresponding  to  the 
desired  index-profile)  it  is  possible  to  establish,  for  para- 
bolic profiles,  an  unstuffing  schedule  which  specifies  the 
time  variation  of  the  concentration  of  the  dopant  compound 
solution  at  the  surface  of  the  rod.  This  is  achieved  by  solving 
| a diffusion  equation  which  describes  the  isothermal  kinetics 

of  the  desorption  of  CsNO^  solution  (inside  the  porous  rod)  for 
cylindrical  geometry. 
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The  unstuffing  (or  desorption)  schedule  developed  by  com- 
puter study  of  the  diffusion  equation  assumes  a constant  diffusion 
coefficient,  D.  This  diffusion  coefficient  characterizes  the 
inter-diffusion  kinetics  in  the  CSNO3-H20  system.  We  have  ob- 
served that  D is  not  a sensitive  function  of  conentration ,C , 

(in  the  temperature  range  of  0 to  100 °C) , but  depends  on  the 
average  pore  size.  It  was  necessary  to  determine  how  the  diffu- 
sion coefficient  depends  on  the  average  pore  size.  If  the  pore 
size  dependence  was  not  strong,  then  a large  pore  size  variation 
could  be  accepted.  On  the  other  hand,  a strong  pore  size  depen- 
dence of  D imposes  strict  demands  on  the  pore  size  distribution 
of  the  porous  glass,  requiring  it  to  be  sharply  peaked. 

A computer  study  was  instituted  to  solve  the  diffusion  equa- 
tions involved  in  the  Molecular  Stuffing  process.  The  goal  of 
this  computer  study  was  to  discover  an  adequate  unstuffing  schedule 
for  the  development  of  parabolic  profiles.  Mathematically,  the 
problem  reduces  to  finding  a proper  time-dependent  boundary  condi- 
tion which  when  combined  with 

(i)  the  diffusion  equation,  and 
(ii)  the  uniform  initial  condition  (t=0  corresponding  to  the 
end  of  the  soaking  step) 

gives  a parabolic  profile  as  a solution.  More  precisely,  let 
C(r,t)  be  the  concentration  of  CSNO3  in  solution  at  some  point  at 
a radial  distance  r from  the  axis  of  the  rod  and  at  a time  t. 

Let  b be  the  radius  of  the  porous  rod,  then  the  initial  condition 
corresponds  to 

c (r, t = 0)  = cQ [0?r* b] 

* 

r • • ..  xr-^ 


(l) 
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and  the  diffusion  equation  (as  expressed  in  cylindrical  geometry)  is 

3c  f 3 2C  1 3d 

3t  “ D l3r2'  r 3r ' (2) 

The  problem  is  to  specify  the  boundary  condition  C(b,t)  = ? so  that 
the  solution  of  (2)  together  with  (1)  and  C(b,t)  will  be  parabolic 
after  a certain  time,  tf. 

This  formulation  approximates  a rod  of  finite  length  as  an  in- 
finite cylinder  which  is  a reasonable  approximation  since  the  ratio 
of  length  to  diameter  for  our  rods  is  greater  than  15.  Another  as- 
sumption employed  in  this  formulation  is  that  the  porous  rod 
soaked  with  CSNO3  solution  can  be  treated  as  a homogeneous  cylinder 
of  CSNO3  solution  where  the  effect  of  porosity  is  accounted  for  in 
the  diffusion  coefficient.  It  is  difficult  to  justify  this  assump- 
tion a priori.  However,  results  from  experiments  have  shown  reason- 
able agreement  with  the  predictions  of  this  formulation  and  support 
this  assumption. 

Carslaw  and  Jaegger  have  derived  a solution  of  Eq.  (2)  cor- 
responding to  the  initial  condition  (1)  and  an  arbitrary  boundary 
condition  C(b,t).  The  solution  is  in  terms  of  an  infinite  series 
in  which  an  integral  of  the  type 

fcf  “n2* 

/ e <p  (A)  d X [where  c(b,t)  = 4>(t)  ] 
o' 

appears  as  a coefficient  of  the  nfc^  term.  It  is  not  clear  whether 
Carslaw  and  Jaegger' s solution  can  be  inverted  analytically  to 
derive  the  proper  <f>(t)  for  parabolic  profiles.  Since  our  primary 
interest  in  not  in  finding  the  particular  boundary  condition  for 
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which  the  solution  will  be  exactly  parabolic,  but  is  to  discover 
an  adequate  boundary  condition  which  can  be  used  in  the  laboratory 
with  relative  ease,  and  for  which  the  solution  can  be  made  as 
parobolic  as  desired  by  adjusting  the  parameters.  We  decided  to 
solve  this  problem  by  using  a computer.  We  have  discovered  an 
unstuffing  schedule  4>(t),  which  gives  nearly  parabolic  profiles. 
These  can  be  made  parabolic  to  any  desired  degree  of  accuracy  by 
proper  choice  of  parameters  shown  in  Fig.  1.  It  shows  C(b,t)  as 
a function  of  unstuffing  time.  Time  t = 0 marks  the  end  of  the 
sorption  step  and  the  beginning  of  the  unstuffing  step.  At  the 
end  of  the  sorption  stage,  the  rod  has  a uniform  concentration  of 
CSNO3  which  we  denote  by  Co*  The  unstuffing  schedule  consists  of 
two  steps: 

(i)  hold : the  surface  concentration  C(b,t)  is  held  fixed  at 
a certain  value  ( Cg)  for  a time,  ti 
(ii)  ramp:  from  time  tj  to  tf  the  surface  concentration  C(b,t) 
is  decreased  linearly  so  that 

t-tj_ 

c (b , t ) = CilWr — i-)]  (3) 

tf-ul 

At  the  end  (t=tf)  of  unstuffing,  the  rod  is  quenched  to  a 
lower  temperature  to  freeze  the  profile. 

The  solution  of  the  diffusion  equation  for  these  boundary 
conditions  can  be  written  in  terms  of  an  infinite  series: 


C(r,tf) 


* 8 Jrt(8  x) 
E 0 n 

n=l  8 3 J (6  ) 
i n 


(4) 
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where 


d> 

y n 


*Vi 


+ eA(Co  - cx)] 


x = r/b 
J0(6n)  = 0 

t — a i 

T . 6 . 


f-Tl 


and  Jq  and  Jf  are  Bessel's  functions  of  zero  and  first  order  re- 
spectively. Notice  that  according  to  Eq.  (4) 

C (b,  tf ) = 0 

00  8 4>n 

and  C (0, tf ) = Z ■ ,R  ' 

n=l  V/n1 

Let  AC  be  the  difference  between  C(0,tf)  and  C(b,tf).  Then 
°°  8 

AC  = E F3J  (6  ) (6) 

n=l  Vl(V 

Notice  that  AC  is  a constant  and  is  determined  by  four  para- 
meters, Cq,  Cf,  tf  and  tf.  Eq.  (4)  can  now  be  rewritten  in  terms  of 
AC  as  follows: 


C (r, t* ) = AC  (l-xz ) + Z 


co  8J0(Bnx) 


[4>n-AC]  (7) 


n=l  Bn’JitBn)  “ 

Eq.  (7)  shows  that  the  resulting  profile  C(r,tf)  is  parabolic 
provided  the  summation  term  is  negligible.  The  higher  order  terms 
in  the  sum  are  indeed  negligible  since  they  contain  in  the  denomi- 
nator a factor  of  8^  and  8n  increases. rapidly  witn  n.  The  first 
three  terms  in  the  sum  can  be  made  to  vanish  by  choosing  the  three 
unknown  parameters  C^,  t^  and  tf  so  that 

AC  = 4'1  = |J>2  = 4>3  (8) 


6. 

Notice  that  is  fixed  by  the  soaking  step.  For  those  values  j 

of  t^  and  tf  which  satisfy  Eq.  (8),  Eq.  (7)  reduces  to 

, 8Jn(S_x) 

C(rft  ) = AC  (1-x2)  + I ~ 3—  [»n-AC] 

f n=4  e3J(0n)  n (9) 

The  assertion  that  the  sum  in  Eq.  (9)  is  negligible  was  verified  j 

on  the  computer  by  actually  summing  up  to  40  terms. 

A typical  simulation  if  Eq.  (9)  by  keeping  40  terms  in  the 

sum  is  shown  in  Fig.  2.  This  profile  corresponds  to  the  following 

values  of  the  parameters: 

Cq  = 0.03,  b = 0.4  cm,  D = 1 x 10-^  cm^/sec. 

= 0.0.4687,  t]_  = 3917  sec,  tf  = 4994  sec. 

C (r,tf ) 

Fig.  2 shows  the  normalized  concentration  — ~ against  normalized 

L0 

radial  distance  x(  = r/b) . The  square  points  represent  the  solution 
of  E<g(9)  and  the  curve  represents  a least  square  parabolic  fit  to 
these  points.  The  deviations,  at  any  radius,  are  on  the  order  of 
10“8  as  shown  in  Fig.  2(b).  Thus  it  is  clear  that  the  unstuffing 
schedule  shown  in  Fig.  1 is  adequate  for  producing  parabolic  pro- 
files. 

As  mentioned  above,  the  unstuffing  schedule  described  in  Fig.l 
is  adequate  for  parabolic  profiling  provided  the  diffusion  coefficient 
D is  constant  and  uniform  throughout  the  rod.  D is  the  coefficient 
of  interdiffusion  of,  for  example,  Cs^O^-I^O  system  inside  the  pores 
of  the  porous  glass.  The  constancy  of  D implies  that  it  is  indepen- 
dent of  concentration  C of  the  solution.  The  uniformity  of  D demands. 
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on  the  other  hand,  the  uniformity  of  the  pore  size  throughout 
the  porous  glass  because  D also  depends  on  the  pore  size.  We  have 
carried  out  experiments  to  study 

(i)  the  concentration  dependence  of  the  diffusion  coefficient 
D,  and 

(ii)  the  uniformity  of  D along  the  axis  of  the  rod. 

The  experimental  technique  to  measure  diffusion  coefficients 
of  solutions  in  the  porous  glass  consists  of  monitoring  the  mass 
gain  as  a function  of  time-  of  a porous  glass  rod  at  a fixed  tem- 
perature in  a CSNO3  solution  of  fixed  concentration. 

A sketch  of  the  apparatus  is  shown  in  Fig.  3.  It  consists 
of  an  aluminum  block  having  the  following  features: 

(i)  a vertical  opening  at  the  center,  having  a diameter  of 

about  1 1/4"  to  hold  a tight  fitting  glass  tube  containing 
a saturated  solution  of  CSNO3 ; 

(ii)  two  smaller  openings  to  place  a thermistor  and  a thermo- 
meter. The  thermistor  is  connected  to  an  external  con- 
troller and  is  used  to  control  the  temperature  of  the 
block; 

(iii)  four  heating  elements  at  the  corners  which  are  connected 
to  a power  supply  via  the  temperature  controller.  The 
temperature  of  the  block  is  held  constant  to  better  than 
0.01°C  and  is  uniform  to  within  0.1°C. 

The  weight  of  the  rod  is  monitored  with  the  help  of  a trans- 
ducer which  converts  mass  variations  into  a voltage  output.  The 
transducer  is  accurate  and  sensitive  to  mass  change  on  the  order 
10~4g  or  larger. 


T' 
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It  can  be  shown  from  diffusion  theory  that  the  time  (t) 
variation  of  the  mass  (M)  of  the  rod  can  be  described  by: 


M (t ) - M (0 ) 

M(oo)  - M ( 0 ) 


(10) 


where  (i)  are  the  roots  of  the  equation  Jq  (Bn)  = 0 (Jq  being 
the  Bessel  function  of  zero  order) . 

(ii)  M(0)  and  M(°°)  are  two  constants  representing  the  mass 
of  the  rod  at  the  start  (t=0)  of  the  sorption  process. 

(iii)  b is  the  radius  of  the  rod. 

Knowing  M(t)  as  a function  of  (t)  it  is  possible  to  evaluate 
D by  a non-linear  least  squares  fit  of  Eq.  (10)  to  the  experimental 
data . 

A typical  measurement  of  M(t)  is  shown  in  Fig.  4.  The  dots 
are  the  experimental  data  points  and  the  curve  is  the  non-linear 
least  squares  fit  of  Eq.  (10)  . It  is  clear  that  the  fit  is  ex- 
cellent. Therefore,  we  believe  that  this  technique  gives  reliable 
estimates  of  the  diffusion  coefficients.  The  diffusion  coefficients 
obtained  in  this  fashion  are  average  diffusion  coefficients  for 
the  rod. 

By  using  solutions  of  different  concentrations,  we  measured 
diffusion  coefficients  as  a function  of  concentration  and  found 
that  D is  relatively  independent  of  concentration  in  the  concentra- 
tion range  of  interest  to  us. 

An  experimental  study  of  diffusion  coefficient  variation 
along  the  length  of  a 12"  rod  was  carried  out  by  dividing  the  rod 
into  six  2"  pieces  and  measuring  the  diffusion  coefficient  in  each 
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piece.  Although  the  cross-section  of  the  rod  was  not  uniform 
along  the  length,  the  results  were  encouraging  in  the  sense  that 
the  variation  of  D between  different  pieces  was  less  than  2%. 

Ill  CONCLUSION 

We  have  demonstrated  by  computer  studies  and  by  experiments 
that  it  is  feasible  to  develop  parabolic  profiles  in  porous  glass 
rods  via  molecular  stuffing.  In  particular  for  parabolic  profiles, 
it  is  important  to  control  the  diffusion  coefficient  inside  the 
pores  to  within  1%.  This  in  turn  requires  a strict  control  of 
glass  composition  and  heat-treatment . We  have  shown  by  experi- 
ment that  such  a control  is  possible. 
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CAPTIONS 


Figure  1. 


Figure  2. 


Unstuffing  schedule  for  parabolic  profiling  by  molecular 
stuffing.  Time  t = 0 marks  the  end  of  soaking  and  the 
beginning  of  unstuffing  (or  desorption)  stage. 

C (r) 

(a)  A plot  of  normalized  concentration  cg  vs. 
normalized  radius  x (=r/b) . The  squares 
represent  the  solutions  of  Eq.  (9)  and  the 
curve  represents  a least  squares  parabolic  fit 
to  the  square  points. 


(b)  The  deviation  between  square  points  and  parabola 
are  shown  as  a function  of  x. 


Figure  3.  A sketch  of  the  apparatus  to  measure  diffusion  coefficient 
of  CsNO-j-^O  system  inside  the  porous  rods. 

Figure  4.  A typical  measurement  of  mass  change  M(t)  as  a function 
of  time  t.  The  dots  represent  experimental  data  and 
the  curve  is  a non-linear  least  squares  fit  of  Eq. (10) . 
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